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ABSTRACT 

Context. About 10 persistently highly absorbed super-giant High-Mass X-ray Binaries (sgHMXB) have been discovered by 
INTEGRAL as bright hard X-ray sources lacking bright X-ray counterparts. Besides IGR J163 18-4848 that features peculiar charac- 
teristics, the other members of this family share many properties with the classical wind-fed sgHMXB systems. 
Aims. Our goal is to understand the specificities of highly absorbed sgHMXB and in particular of the companion stellar wind, thought 
to be responsible for the strong absorption. 

Methods. We have monitored IGR J17252 - 3616, a highly absorbed system featuring eclipses, with XMM-Newton to study the vari- 
ability of the column density and of the Fe Ka emission line along the orbit and during the eclipses. We also built a 3D model of the 
structure of the stellar wind to reproduce the observed variability. 

Results. We first derived a refined orbital solution for this system built from INTEGRAL, RXTE and XMM-Newton data. The XMM- 
Newton monitoring campaign revealed significant variation of intrinsic absorbing column density along the orbit and of the Fe Ko- line 
equivalent width around the eclipses. The origin of the soft X-ray absorption is modeled with an dense and extended hydrodynamical 
tail, trailing the neutron star. This structure extends along most of the orbit, indicating that the stellar wind is strongly disrupted by 
the neutron star. The variability of the absorbing column density suggests that the wind velocity is smaller (v^ « 400 km/s) than 
observed in classical systems. This can also explain the much stronger density perturbation inferred from the observations. Most of 
the Fe Ko- emission is generated in the most inner region of the hydrodynamical tail. This region, that extends over a few accretion 
radii, is ionized and does not contribute to the soft X-ray absorption. 

Conclusions. We have built a qualitative model of the stellar wind of IGR J17252 - 3616 that can represent the observations and 
suggest that highly absorbed systems have a lower wind velocity than classical sgHMXB. This proposal could be tested with de- 
tailed numerical simulations and high-resolution infrared/optical observations. If confirmed, it may turn out that half of the persistent 
sgHMXB have low stellar wind speeds. 

Key words. X-Rays:binaries - Stars: pulsars : individuals: IGR J17252 - 3616= EXO 1722-363 - 

1 . Introduction Highly absorbe d sgHMXB were discovered by INTEGRAL 

dWalteretaUl2004 and are characterized by strong and persis- 

High mass X-ray binaries (HMXB) consist of a neutron star or tent soft X-ray absorption (N H > 10 23 cm 4 ). When detected, 

a black-hole fueled by the accretion of the wind of an early type these systems have short orbital periods and long spin periods 

stellar companion. Their X-ray emission, a measure of the ac- (IWalter et al]l2006l) . They correspo nds to the cat egory of wind- 

cretion rate, shows a variety of transient to persistent patterns, fed accretors in the Corbet diagram dCorbetlll986h . 

Outbursts are observed on timescales from seconds to months lr ^ n T1 „_ c _ , u . . , * . rv^nr u a 

. . . . , „ . ,_ 4 _ . . IGR J17252 - 3616 has been detected by IS GRI on board 

and dynamical ranges varying by factors of 1 C . The majorit y mTEG RAL Qn Februarv 9 , 2004 a other hard X -ray 

of the known HMXB are Be/X-ray binaries <Lm et alj l2QQfl>, sources Walter et al.ll2004l [2006). The source has first been de- 

i3 , with Be stellar companions. These systems are transient, featur- . . , , r vr\c a-t ,rvr>. n-n ^hml\ i ^ v 

C3 . , . . . . . . . . . , , J , tected by EXOSAT (EXO 1722-3616) as a weak soft X-ray 

ing bright outburst with typical duration on the order of several / , . lri0/1 irrr — i ii inoot , t inc? ^- 

weeks Icharles &C oe 2006). A second class of HMXBs har- source back in 1 984 ^Warwick et al- EHi). In 1987 Gtnga per- 

: 131 . . ™„„, , „ ,. , formed a pointed observations and revealed a highly variable 

bour OB supergiants companions (sgHMXBs) feeding the com- X -ray source, X17 22-363, with a pulsations of ~ 413.9 sec 

pact object through strong, radiatively driven, stellar winds or — ' ' , | ln0 n k t- .u ^ ■ u »■ i j 

£. , , . a m . ,.,„ r i nl , , , , (Tawara et al. 1989). Further Ginga observations revealed the or- 

Roche lobe overflow. Thanks to INTEGRAL, the number known rr - ; . „ ° . 

T „„,_ ■ , , , • , , 1 jTTT - ; n bital period of 9-10 days and a mass of the companion star of ~ 

sgHMXB systems tripled during the last decade « Walter et alj ^ Mq J T akeuchiet alJll990h . Both papers concluded that the 

system was a high mass X-ray binary (HMXB). 

INTEGRAL and XMM observations of IGR J17252 - 



Send offprint requests to: A. Manousakis 36 fallowed to identify the infrared counterpart of the system 

* Based on observations obtained with XMM-Newton and and t0 accurately measure the absorbing column density and 

INTEGRAL, two ESA science mission with instruments, data refine the spin period of the system. Thanks to the eclipses, 

centers, and contributions directly funded by ESA Member States, an ac curate orbital period cou ld be derived from INTEGRAL 

NASA, and Russia. data dZurita Heras et al 1 120061) . Further RXTE observations re- 
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vealed a highly inclined system (; > 61°) w ith a companion 
star o f M, < 20 M G and R, ~ 20 - 40 R Q (iThompson et al.l 
120071) . Re cent VLT observa t ions provided the co mpanion spec- 
tral type dChatv et al.ll2008t iMason et alJl2009al) and radial ve- 
locities measurements (IMason et alj|2009bl) . Its spectral energy 
distribution can be characterize d by a temperature of T* ~ 30 
kK and a reddening of Ay ~ 20 (iRahoui et al.ll2008l) . 

In this paper, we report on a monitoring campaign of IGR 
J 17252 - 3616 performed with XMM-Newton along the orbit 
in order to estimate the structure of the stellar wind and of the 
absorbing material in the system. We describe the data and their 
analysis in Section 2, a refined orbital solution in section 3 and 
present the evolution of the X-ray spectral shape along the orbit 
in section 4. In section 5, we present and discuss a 3D model of 
the structure of the stellar wind that can reproduce the observa- 
tions and conclude in section 6. 



2. Data reduction and analysis 

2.1. XMM-Newton 

Pointed observations of IGR J17252 - 3616 were per- 
formed between A ugust and October 2006 with XMM-Newton 
dJansen et al. I l200l . We scheduled 9 observations to cover the 
orbital phases : 0.01, 0.03, 0.08, 0.15, 0.27, 0.40, 0.65, 0.79, and 
0.91. In addition, we used one observation of 2004 with a phase 
of 0.37. The observations are summarized in tableQ] 

The Science Analysis Software (XMM-SAS) version 9.O.C0 
was used to produ ce event lists for the EPIC-pn instrument 
dStruderetaljboOll) by running epchain. Barycentric correc- 
tion and good time intervals (GTI) were applied. Photon pile- 
up and/or out-of-time event were not identified among the data. 
High level products (i.e., spectra and lightcurves) were produced 
using evselectQ. Spectra and lightcurves were built by collect- 
ing double and single events in the energy range 0.2 - 10 keV. 
The lightcurves were built using 5 sec time bins. The spectra 
were re-binned to obtain 25 counts/bin for low count rate obser- 
vations and 100 counts/bin for high count rate observation. 

In one dataset (ObsID 0405640701, ~ 19 ksec), the count 
rate above 10 keV has a very peculiar behavior, increasing 
monotonically with time. As this does not affect the background 
subtracted source lightcurve significantly we kept the entire data 
for pulse arrival times determination. Standard GTI was used for 
spectral analysis. 

2.2. INTEGRAL 

We reanalysed the complete hard X-r ay light-curve of I GR 
J17252 - 3616 obt ained with ISGRI dLebrun et alJ 120031) on 
board INTEGRAL dWinkler et alj|2003l) . We extracted the 22- 
40 k eV light-curve using the HEAVENS interface dWalter 
120 101) provided by the INTEGRAL Science Data Centre^ 
dCourvoisier et alj |2003). The lightcurve includes all public data 
available on IGR J17252 - 3616 from 2003-01-29 06:00:00 
to 2009-04-08 00:28:48 UTC. The effective exposure time on 
source is ~ 3.6 Msec. 



1 http://xmm.esac.esa.int/sas/ 

2 http://xmm.esac.esa.int/sas/current/howtousesas.shtml 
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3. Timing analysis and orbit determination 



3.1. Orbital Period from INTEGRAL 



We have used the Lomb-Scarg le (lPress&Rvbickilll989h tech- 
nique to determine the orbital period from the INTEGRAL light- 
curve and obtained P ,-£=9.742 + 0.001 days. Figure Q] shows 
(upper panel) the Lomb-Scargle power around the orbital period 
(dashed line). This period was used to refine the orbital solu- 
tion (section [33T l. The lower panel of figure [T] shows the light- 
curve folded with the newly derived orbital period. The eclipse 
is clearly detected with the count rate dropping to zero. 
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Fig. 1. Top: Lomb-Scargle periodogram obtained from the 
INTEGRAL 22-40 keV lightcurve. Bottom: INTEGRAL 
lightcurve folded with a period of 9.742 days. 



http://www.isdc.unige.ch/heavens 
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Table 1. XMM-Newton observation log listing the revolution number, the start time, the effective exposure, the source counts, and 
the phase, calculated at the middle of each observation using the orbital solution obtained with fixed orbital period (table[3j. 



ObsID 


# 


Revolution 


Start time (UT) 


Effective 


Source counts 


Phase 










exposure (ks) 




(+ 0.01) 


0405640701 


I 


1231 


7006-08-79X0^-07 -58 


19.2 


5.4xl0 2 


0.03 


0405640301 


2 


1232 


2006-08-31T16:37:44 


4.1 


1.4xl0 4 


0.27 


0405640401 


3 


1234 


2006-09-04T06:35:33 


5.6 


7.9xl0 3 


0.65 


0405640501 


4 


1235 


2006-09-06T19:33:10 


5.4 


9.5xl0 2 


0.91 


0405640601 


5 


1236 


2006-09-08T10:03:38 


7.9 


3.8xl0 2 


0.08 


0405640701 


6 


1239 


2006-09-15T07:23:45 


2.8 


1.7xl0 3 


0.79 


0405640801 


9 


1247 


2006-10-01T03:24:26 


9.4 


2.3xl0 4 


0.40 


0405640901 


8 


1246 


2006-09-28T14:36:53 


11.3 


1.9xl0 4 


0.15 


0405641001 


7 


1245 


2006-09-27T07:27:58 


9.4 


6.5xl0 2 


0.01 


0206380401 


10 


785 


2004-03-21T13:23:09 


8.6 


4.7xl0 4 


0.37 



Table 2. Pulse Arrival Times and derived pulse period. Last col- 
umn gives the Lomb-Scargle periods. We couldn't get a good 
S/N LS diagram for observation 0405640301. 



ObsID PAT(HJD) Pulse period (s) L-S period (s) 



(±0.00011 ) 



0405640301 


53978.74948 


414.3 : 


tO.l 






0405640401 


53982.32329 


414.0 : 


tO.l 


417.4 : 


t0.6 


0405640701 


53993.45320 


414.2 : 


t 0.2 


414.6 : 


tO.l 


0405640801 


54009.24191 


414.2 : 


tO.l 


414.2 : 


t0.4 


0405640901 


54006.75737 


413.8 : 


tO.l 


413.8 : 


t0.3 



3.2. Pulse Arrival Times 

Pulse arrival times (PATs, hereafter) have been obtained from the 
broad-band 0.2-10 keV lightcurves obtained by XMM-Newton. 
Close to the eclipse (0 = 0.03, 0.08, 0.91, 0.01), when the com- 
pact object is behind the massive star, pulses could not be de- 
tected. We did not extract PATs for the observation of 2004. 

In order to determine the PATs we used a pulse profile tem- 
plate. This template is derived by folding the lightcurve from 
observation 0405640801 with a p eriod of 414.2 sec, ob tained 
using the Lomb-Scargle technique dPress & Rvbic ki 1989). This 
observation was selected because the source was very bright for 
a long and almost uninterrupted exposure. 

A sequence of pulse profile template was fit to each individ- 
ual lightcurve. This sequence is characterized by: (i) the time of a 
pulse at the middle of the observation, (ii) the pulse period, and 
(iii) the amplitudes of each pulse. This assumes that the pulse 
period is reasonably constant during each observation. 

The statistical errors on the pulse time and period are typi- 
cally 0.01 sec and 0.1 sec, respectively. Table [2] lists the pulse 
times and the spin period obtained at the middle of each obser- 
vation. The spin periods were always consistent with the Lomb- 
Scargle period of each run. 

The PAT accuracy is however limited by the systematic error 
related to the assumed pulse template. Using a different pulse 
profile template (derived from observation 0405640901) pro- 
duces PAT with offset of 8 to 12 sec from the values listed in 
table|2] We adopted a systematic error of 10 sec. 

3.3. Orbital Solution 

We have derived the orbi tal solution using t he PAT s of the RXTE 
observation obtained by iThompson et alj ( 120071 Epoch 3), and 
the PATs derived above from XMM data. 



The orbital solution is obtained by comparing the observed 
pulse arrival delays (t„ -to- nPp - \n 2 PoP) to th e expected ones 
ff. r sin/ cos[27r(f„-r9o)/fVz>] dLevine et al.l2004l) . The orbital pa- 
rameters (the orbital period, P,„b\ the projected semi-major axis, 
a x smi; the reference time corresponding to mid-eclipse, Tgo) are 
assumed to be constant. To account for pulse evolution, two 
set of pulse parameters (spin period at time to, Po', spin period 
derivative, P) have been used for the RXTE and XMM cam- 
paigns. The pulse number n is given by the nearest integer to 
n = (f„ - f )/(P + 0.5P(f„ - to)). 

We have performed a combined fit for the RXTE and XMM 
observations. We derived the orbital solutions by (i) allowing all 
the parameters to vary freely and (ii) by fixing the orbital period 
to the value derived from INTEGRAL data. The resulting pa- 
rameters are listed in table [3] We could also obtained an upper 
limit (90%) of e < 0.15 on the eccentricity. The latter has been 
achieved by adding the fir st-order term in a Taylor series ex- 
pansion in the eccentricity dLevine et al 1 12004 . The RXTE and 
XMM orbital solutions are comparable and the resulting param- 
eters are consistent within the errors. 

The folded lightcurve obtained by INTEGRAL derived or- 
bital period (fig. [1] lower panel) results in a pulse fraction 
~100%. For the rest of the analysis, we used the orbital solu- 
tion obtained with fixed P ,i, (table[3j. 

Figure [2] shows the resulting pulse arrival times delays (fixed 
Porb) for both RXTE and XMM data together with the best-fit 
orbital solution. 



4. Spectral analysis 

The spectral analysis was performed using the XSPEC0 pack- 
age version 11.3.2ag (lArnau dl 119961) . In order to use the xl 
statistics, we grouped the data to have at least 25 (faint 
spectra) and up to 100 (bright spectra) counts per bin. We 
initially fitted the observed spectra using a phenomenolog- 
ical model made of an intrinsically absorbed cutoff power 
law, a blackbody soft excess, and a gaussian Fe Ka line 
(wabs*(bb+gauss+vphabs*cutoff)). The centroid of the 
iron line is compatible with E c ~ 6.40 + 0.03 keV at all epochs. 

The spectra are always strongly absorbed below ~ 3 keV and 
feature an iron K-edge at ~ 7.2 + 0.2 keV. We first fit each spec- 
trum with all parameters free, e xcepting the galactic absor ption, 
fixed to N H = 1.5 x 10 22 cirT 2 (iDickev & Lockmanll 19901) . 



4 http://xspec.gsfc.nasa.gov 
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Fig. 2. Del ays (top) and residuals (bottom) derived from the fixed P or b orbital solution compared to the data. Left: RXTE data are 
taken from lThompson et al.l (120071) . Right: XMM - Newton data from this work. 



Table 3. Orbital solution for IGR J17252 - 3616. The errors 
have been calculated at 90% confidence level. The errors on 
the arrival times t are 0.00001 and 0.0001 days for RXTE 
and XMM, respectively . The last column shows the result from 
iThompson et al.1 (120071) for comparison. 





Units 


Free P Brb 


Fixed P orb 


RXTE Epoch 3 


,XTE 
pXTE 
pXTE 


HJD 

s 

fIS s~' 


53761.73144 
413.889± 0.004 
-0.010±0.002 


53761.73142 

413.889±0.005 

-0.010±0.003 


53761.73126 
413.894± 0.002 
-0.0106±0.0001 


t XMM 
pXMM 
pXMM 


HJD 

s 

fJS s _1 


53978.7494 
413.86 ±0.04 
0.98± 0.04 


53978.7495 
413.84± 0.04 
1.01 ±0.03 




o^-sim 

Porb 


lt-s 
d 

HJD 


102±8 
9.76 ± 0.02 
53761.62 ±0.1 


101±2 

9.742 (fixed) 
53761.69+0.1 


101±4 
9.78±0.04 
53761.60±0.09 


xl 




1.5 (30) 


1.65 (31) 


1.45 (28) 



2 0.01 =• 




Energy (keV) 

Fig. 3. Top: Folded model and data for 0=0.65 (black) and 
0=0.91 (red). Bottom: The corresponding residuals using the 
model described in the text with the parameters from table [4] 

Some parameters (photon index, cutoff energy, blackbody 
temperature, absorber Fe metallicity) did not vary (within the 
90% errors) among the observations and were fixed to their av- 
erage values of Ec = 8.2 keV, F = 0.02, and kTBB=0.5 keV, Z=l Z . The Fe line energy is fixed to E=6.4 keV, with a nar- 
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row width. Two representative spectra of our observations are 
displayed in figure [3] 

The intrinsic absorbing column density and the normaliza- 
tion of each component could vary freely. The best fit parameters 
are listed in table|4] Figure|4]shows the variation of all these pa- 
rameters and of the unabsorbed 2-10 keV flux and Fe Ka equiv- 
alent width (EW). The unabsorbed Fe Ka EW was calculated 
by setting the intrinsic absorbing column density to zero. All the 
spectral fits were good resulting in xl ~ 0.8 - 1 .4, apart for 
the 0=0.03 observation providing a poor fit. During the eclipse, 
some parameters are poorly constrained. 

The unabsorbed continuum flux (fig. 3rf) is of the order of 
~ 5 x 10~" erg cirT 2 s outside of the eclipse. Variations are 
observed between (2 - 10) xlO 11 erg cirT 2 s _1 and they can 
be interpreted as variation of the accretion rate (M). During the 
eclipse the continuum flux drops by a factor of ~ 200 and re- 



mains stable at the level of F 



2-Wkev 
unabs 



(7 + 1) x 10~ 13 erg cm 



,-2 



The intrinsic absorbing column density (fig. |4j-b) is persis- 
tently high (<; 10 23 cm -2 ). Significant variations are detected 
for <p = 0.2 - 0.4 and close to the eclipse, reaching values of 
N H ~ 9 x 10 23 crrr 2 . 

The normalization of the blackbody component (fig-H-c) 
does not show any variability, although the low energy part of the 
spectrum is poorly constrained. The normalization of the black- 
body component is compatible with ~ 10 33 erg s _1 assuming 
a distance of 8 kpc. The very high intrinsic absorbing column 
density rules out the neutron star as the origin of the soft excess. 

The flux and EW of the Fe Ka line are displayed on fig- 
ure [4j-d and]4j-e, respectively. Both components show significant 
variations indicating that the region emitting the line is at least 
partially obscured by the mass-donor star. 

As the X-ray continuum illuminating the gas emitting the Fe 
fluorescent line cannot be measured during the eclipse we calcu- 
lated a corrected Fe Ka EW by assuming a constant continuum 
flux of 1.8 x 10~ 3 ph keV _1 cnT 2 s" 1 (fig.0. 

5. Discussion 

5. 1 . Constraining the physical parameters of the system. 

In the previous sections we have derived an orbital solution (ta- 
ble [3} yielding a mass function / = 4n 2 (a x sin i) 3 /GP 2 rb = 
M Bsin i il(\ + q) 2 « 11.7 + 0.7 M Q . This is consistent with 
a high mass x-ray binary system. 

Adopting an inclination ; = 90° we infer a mass Mob ~ 14 
M for the dono r star. Radial velocit y observation showed q = 
M x I M b ~ 0. 1 dMason et al.ll2009b(). Using an upper limit on 
the mass of ~ 2OM ( Thompson et al. 2003 constrain tne m ~ 
clination of the orbit i > 70°. iMason et alJ (l2009bl) estimated 
an inclination ; as 75° - 90°. With this value the mass of the 
donor star is constrained in the range Mob ~ 14-17 M . 
The mass ratio and the donor mass imply a neutron star mass 
Mns = 1.4—1.7 M . Both the masses of the donor star and of 
the comp act object are roughly similar (within a factor of 2) to 
VelaX-1 (lOuaintrell et"alll2003l:lvan Kerkwiik et aT]ll995l) . 

The separation of the system could be derived from the dura - 
tion of the eclipse and the inclination (Joss &Rappaportlll984l) . 
For the eclipse duration of A0 * 0.18 + 0.02 (i.e. 1.75 days) we 
can estimate the separation a x ss 1.7 — 1.8/?*. The Roche lobe 
(/?/.) of the system sc ales from Ri - 0.99 - 1 . 06 R* assuming a 
synchronous rotation dJoss & RappaporJll984l) . This means that 
the system is very close to fill its Roche lobe and form an accre- 
tion disk, although no significant spin-up has been observed. 



Using our VLT observations, Mas on et al ] (|2009a) per- 
formed near-IR spectroscopy of IGR J 17252 - 3616 and con- 
cluded that the donor star is a B0-B5 I to B0-B1 la with effective 
temperature T e ff = 22-28 kT and a stellar radius of/?* = 22-36 
R . Based on a spectroscopic determination, we can translate it 
into absolute numbers for the separation (a x « 37 - 64 R ) and 
the Roche lobe radius (/?£ « 22 - 38 R ). 

The unabsorbed 2-10 keV source flux is in the range (0.2 - 
1.3) x 10~ 10 erg s _1 cirT 2 . Adopting a mean value o f (0.8 + 0.3) x 
10 10 er g s _1 cm~ 2 and assuming a distance of 8 kpc dMason et al.l 
l2009al) the inferred 2-10 keV unabsorbed luminosity is Lx ~ 
10 36 erg s . Assuming accretion as the source of energy (Lx = 
eMc 2 , we can estimate a mass accretion rate of M ~ 10~ 9 M 
yr" 1 similar to that of Vela X-l dFiirst et al.ll2010h . 

During the eclipse, the X-ray luminosity drops by a factor of 
~ 200 resulting a Lx ~ 5 x 10 33 erg s _1 . As OB st ars are emitting 
in X-r ays with a luminosity Lx ~ 10 3I ~ 32 erg s _1 dGiidel & Nazel 
120091) . this e mission is pro bably dominated by scattering in the 
stel lar wind (lHaberlll 19911) . 

iHickox et al.1 (|2004) have discussed the origin of soft X-ray 
excesses in many types of accreting pulsars. It is likely, that the 
fairly constant soft X-ray excess o bserved in IGR J 17 252 - 3616 
is emitted by recombination lines dSchulz et al.l2002|) in a regio n 
of the wind larger than the stellar radius dWatanabe et al.ll2Q06b . 

For a spherically symmetric stellar wind, one would expect 
to have smooth and predictable variations on the absorbing col- 
umn density along the orbit. In particular observations on the 
same line of sight or symmetric when compared to the eclipse 
must result in different, respectively identical, column densities. 
Our observing strategy resulted in Nu(<p -0.15) ~ Nui<p = 0.37), 
N H ((f> = -0.35)> N H (0 = 0.37) and N H (</> = -0.35)* N H (<f> = - 
0.21), ruling out a spherical geometry. 



5.2. Stellar wind structure 

We have constructed a 3D model of the OB supergiant stellar 
wind in order to simulate the variability of the intrinsic column 
density (Nh) and of the Fe Ka line equivalent width (EW), along 
the orbit, and to compare them with the observation s. We have 
assumed a distance D * 8 kpc dMason et al ] l2009al) . a circular 
orbit (e = 0) and an edge-on geometry (i = 90°). 



Variability of N H along the orbit 

We first investigate the behavior of the intrinsic column density, 
Nh, as a function of phase, to reveal the structure of the wind 
along the orbit. We approximated the wind structure with two 
components, the unperturbed wind (p, v ,„ ( /) and a tail-like hydro- 
dynamic perturbation (p, a u) related to the presence of the neutron 
star. Ta il-like structures are predicted by hydrodynamical simu- 
lations dBlondin et alj|1990t) but produce N H , up to ~10 22 cm 2 , 
too small to account for the variability observed in IGR J17252- 
3616. 

The unperturbed ste llar wind has been modeled by assuming 
a standard wind profile dCastor et al.ll 19751) . 

/?*y 



/ /?,V 
v(r) = fool 1 - —J 



where v{r) is the wind velocity at distance r from the stellar cen- 
ter, Voc is the terminal velocity of the wind, and f3 is a parame- 
ter describing the wind gradient. The conservation of mass pro- 
vides the radial density distribution of the stellar wind. The un- 
perturbed stellar wind is a good approximation within the orbit 
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Fig. 4. Spectral variability along the orbit, a) The cut-off power-law normalization at 1 keV, b) The intrinsic hydrogen column 
density, c) The soft excess blackbody luminosity, assuming a distance of 8 kpc . d) The Iron line flux, e) The iron line equivalent 
width calculated for the unabsorbed continuum, and f) the unabsorbed flux. Dashed vertical lines indicates the eclipse boundary. 
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Table 4. Spectral analysis results. All the free parameters are listed below together with unabsorbed 2-10 keV flux and;^ 2 . All the 
uncertainties have been calculated at 90% confidence level. For more details see the text. 



Phase 


N« 


unabs EW (Fe) 


Fe line flux 


Blackbody Luminosity 


jcutoftpl (1 kf;V) 


punabs 
r 2-10keV 


xl (dof) 


(± 0.01) 


i r.22 -2 

10- z cm - 


keV 


ir\-5 i -1-1 

10 ph cm 1 s ' 


10 33 ergs s~' @ 8 kpc 


10~ 3 ph keV -1 cm~ 2 s" 1 


ir>-!0 -1 -2 

10 erg s cm 




0.01 


i 7 +14 

11 -9 


0.50±0.4 


().34±0.21 


< 0.88 


015 +,mw 


0.0054 ± 0.003 


1.34 (23) 


U.UJ 


26 -6 


U.o 1 +U.4 


n fiA+012 

0.66_ 015 


< Yj.jd 


n n 1 8+0.009 
0.018 

023 +i ">' 2 


U. 000(0)+ U.UUj 


Z.o \ Y 1 ) 


0.08 




1.1 ±0.6 


1.1 ± 0.3 


< 0.57 


0.008(8) ± 0.003 


1.06 (11) 


0.15 


15.7 ±0.5 


0.04 ±0.01 


2 6 +0 - 8 


1.6± 1.4 


1.44 ±0.03 


0.48 ± 0.02 


0.90 (172) 


0.27 


Q 4+0.4 


0.14±0.02 


14± 02 


4.4+0.3 

3 4+°' :2 


2.29± 0.05 


0.77 ± 0.02 


1.05 (128) 


0.37 


12.1+0.3 


0.09±0.006 


16+J 


rt o£+0.05 
3.8D_ 006 


1.29 ±0.04 


1.09 (396) 


0.40 


19.1+0.5 


0.11+0.01 


11±1 


2.8 1+° 06 


2.26±0.05 


0.76 ± 0.02 


1.45 (207) 


0.65 


24+1 


0.11 ±0.02 


7 6+ - 9 
-15 


1 S7°l' 


i cr+0.06 
-0.05 


0.52 ± 0.02 


0.99 (73) 


0.79 


23±3 


0.09±0.03 


-0.7 


< 1.48 


0.60 ±0.06 


0.20 ± 0.04 


0.97 (31) 


0.91 


on+ll 
0, -14 


0.12±0.03 


•3 9+O.5 
J - z -0.7 


1 41 +0 1 

1- -0.05 


6 +0 - 2 


0.20 ±0.07 


0.82 (33) 



of the neutron star. Hydrodynamical simulations 



199011 199 11: iBlondinlll 994 fBlondin & Woolll995t 



Blondin et al 



Mauche et al 



2008) of HMXB have shown that the wind can be highly dis- 



rupted by the neutron star beyond the orbit. 

To estimate the terminal velocity of the unperturbed wind, 
we have studied the Nh variability using three different sets of 
parameters (Pig. |5j. The mass-loss rate and terminal velocity are 
constrained by the data to be in the range M H ,/voa ~ (0.7 - 2) x 
10~ 16 M / km; (f3 has a very limited impact on the results, we 
used 0.7). 

The fraction of the wind captured by the neutron star can be 
estimated from the accretion radius r acc = 2GM x liir )rb + v 2 ) ~ 
2 ■ 10 11 cm (where v or t, = 250 km/s is the orbital velocity) as 
/ ~ nrf KC /4-nR 2 rh ~ 7.5 • 10™ 4 . The mass-loss rate is therefore 
M w ~ f~ l M « 1.5 • 10~ 6 M /yr and the terminal velocity of the 
wind constrained in the interval u m ~ 250 - 600 km/h. 

In our simulation we have adopted a terminal velocity 
foo=400 km/sec, a stellar radius R„-29 R e , a wind gradient 
j8=0.7, and a mass loss rate M*= 1.35xl0 -6 M yr~'. A sum- 
mary of the assumed, observed, and inferred parameters of the 
model is listed in table [5] The variability of Nh along the or- 
bit indicates that the unperturbed wind is adequate for phases 
^0 - 0.35. For <p Z 0.35, N H is increasing by ~ 2 x 10 23 
cm -2 . This indicates that a high density tail-like structure lies 
on one side of the orbit, trailing the neutron star. The tail-like 
component is still present up to ^ 0.8. 

We have assumed that, the tail-like structure is created very 
close to the NS and is opening with distance from the neutron 
star. The density of the material inside the 'tail' is decreasing 
with distance for mass conservation. Its distribution follows a 
'horn' -like shape with a circular section. We have adjusted the 
density of the tail-like structure to match the observations. The 
density distribution p wm d + ptaii is displayed in figure [6] The su- 
pergiant is located at the center (black disk). The tail-like struc- 
ture covers about half of the orbit. 

Figure|5]displays the simulated Invariability from the above 
density distribution together with the observed data points. The 
data and the model shows that the tail-like perturbation is essen- 
tial to understand the observed variations of the Nh- 



Variability of the Fe Ka line along the orbit 

Assuming that the intrinsic X-ray flux is not affected by the 
eclipse, the Fe Ka equivalent width drops by a factor ~ 10 during 
the eclipse in an orbital phase interval of ~ 0.1. This indicates 
that the radius of the region emitting Fe Ka is smaller than half 




0.2 



0.4 0.6 

Phase 



0.8 



Fig. 5. Simulated Nh variations plotted together with the data. 
We illustrate three different wind configurations, for M/Voa = 
0.7 (cyan), 1 (green), and 2 (red) xl0~ 16 M Q /km, keeping all 
the other parameters fixed. The solid black line shows the total 
Nh consisting of the unperturbed stellar wind (green line) and 
of the tail-like extended component. The observations during the 
eclipse have been omitted. 



of the stellar radius (< 10 12 cm) and much more compact than 
the tail structure responsible for the Nh variability profile. 

Vela X- 1 shows a similar beha vior that was interpr eted with 
an emi tting region of th e size of dOhashi et alj|1984h . or even 
within dEndo et al.| [2002). the accretion radius. 

Outside of the eclipse, the equiva lent w i dth of the Fe Ka line 
is of the order of 100 eV. Following iMatB J2002h and assuming 
a spherical transmission geometry, this corresponds to a column 
density of Nh ~ 2 • 10 23 crrT 2 . As this additional absorption is 
not observed, the region emitting Fe Ka is partially ionized. 

A ionization parameter £ = L/nR 2 in the range 10-300 is 
required to fully ionize light elements contributing to the soft X- 
ray absorption and to keep an Fe Ka line at the en ergy of 6.4 keV 
dKallman et al.ll2004lKaliman & McCravlll982h . The density of 
the Fe Ka emitting region is therefore ~ f 'Ri2~ 2 ■ 10 12 cm , 
where R12 is the distance from the neutron star in units of 10 12 
cm. As N H = nR ~ 2 x 10 23 cirT 2 we have £ ~ 5/Rn. A dense 
cocoon is therefore needed around the neutron star with a size 
0.5 > R12 > 0.02 and a density 10" cirT 3 < n < 3 ■ 10 12 cnr 3 . 
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Fig. 6. Number density distribution in the plane of the orbit in- 
cluding a smooth wind and a tail-like perturbation. The black Fig. 7. Corrected unabsorbed Fe Ka line equivalent width along 
disk at the center represents the supergiant star. the orbit. The two curves represent two different configurations 

in the presence of the hydrodynamical structure. Black curve: 
the central cocoon is present. Green curve: without the central 
cocoon. See the text for more details. 



To match the observed column density variability, the radius 
and the density of the inner region of the tail structure were set 
to 4 • 10 11 cm and 3 • 10 11 cm -3 respectively. It is therefore very 
likely that the dense cocoon corresponds to the inner and ionized 
region of the hydrodynamical tail. 

We have thus added this partially ionized cocoon in our sim- 
ulations, using a density of 3 • 10 11 cm -3 within a radius of 
~ 6 ■ 10 11 cm ~ 3R acc . The Fe Ka emissivity map (Fig. [8} is 
built applying an illuminating radiation field (~ 1/r 2 ) to the den- 
sity distribution. 

Figure [7] displays the resulting simulated profile of the Fe 
Ka equivalent width together with the observed data. The green 
curve shows the variations of the Fe Ka equivalent width ex- 
pected from the wind density profile excluding the central co- 
coon, which obviously could not reproduce the data. The black 
curve accounts for the dense central cocoon. The exact profile of 
the eclipse is related to the size and density profile of the cocoon. 
No effort has been made to obtain an exact match to the data. 

The ionized cocoon is expected to produce an iron K-edge 
at ~ 7.8 keV. For a column densit y of N H ~ 2 x 10 23 c m" 1 , the 
optical depth of this edge r ~ 0.2 jKallman et al. 2004) remains 
difficult to detect. Even our observation at an orbital phase 0.14 
(the best candidate for the detection of the ionised edge) does 
not have enough signal. 

The mass of the tail-like structure M ta a ~ 1O~ 8 M can be 
accumulated in t ta u = M lai i/M, where M = {nr 2 ej .jlAnR 2 orb )M w , 
where r e ff is an effective radius for the funneling of the wind in 
the tail. For a tail accumulation time scale (t ta n) comparable to 
the orbital period of (~ 10 days), this effective radius becomes 

12 Race 

The orientation of the tail-like structure depends on the wind 
and orbital velocity. The angle between the wind velocity and the 
orbital velocity is given by tan(a) = v orD / v. The tail obtained in 
our simulations is tilted by a ~ 80°. This would correspond to 
v ~ 0.2 v or b ~ 50 km/s which is lower than v{R,„b) ~ 250 km/s 
because of the ionization of the stellar wind in the vicinity of the 
neutron star. 




Fig. 8. Integrated Fe Ka emissivity centerred on the neutron star 
at phase <p = 0.5 (relative units). The smooth circular halo shows 
the rim of the supergiant star. The tail structure can be observed 
on the right. 



6. Conclusions 

We have presented the analysis of an observing campaign per- 
formed with XMM - Newton on the persistently absorbed sgH- 
MXB IGR J17252 - 3616. Nine observations have been per- 
formed over about four weeks, distributed on various orbital 
phases. Three of them were scheduled during the eclipse of the 
neutron star by the companion star. 

We first refined the orbital solution, using in addition archival 
INTEGRAL and RXTE data and found an orbital period of 
9.742 d and a projected orbital radius of 101 + 2 lt-s. The pul- 
sar spin period varies between 414.3 and 413.8 sec during the 
observing campaign. 

The X-ray spectrum (0.2 - 10 keV), which varies along the 
orbit, was successfully fit using an absorbed cut-off power law 
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Table 5. Summary of the wind model used for IGR J 17252 - 
3616. 



Parameter 


Value 


Reference 


Observed 


q=M x /M 0B 


0.1 


Mason et al. (2009b) 


M 0B (M ) 


15 


Takeuchi et al. (1990) 


a(R„) 


1.75 


This work. 


M (Mo/yr) 


10" 9 


This work. 


Assumed 


R. (R Q ) 


29 




M x (M Q ) 


1.5 




i (deg) 


90 




P 


0.7 




Inferred 


Uoo (km/s) 


400 




(Mo/yr) 


1.35 x 10~ 6 





continuum, a soft excess and a gaussian emission line. The soft 
excess, modeled with a black body, remained constant. 

The continuum component varies in intensity (a measure of 
the instantaneous accretion rate) but features a constant spectral 
shape, as usually observed in accreting pulsars. 

The absorbing column density and the Fe Ka emission line 
show remarkable variations. The column density , always above 
10 23 cm -2 , increases towards 10 24 close to the eclipse, as ex- 
pected for a spherically symmetric wind. The wind velocity is 
unusually small with = 400 km/s. An additional excess of 
absorption of 2 ■ 10 23 cm~ 2 is observed for orbital phases <p > 0.3 
revealing an hydrodynamical tail, trailing the neutron star. 

During the eclipse, the equivalent width of the Fe Ka line 
drops by a factor > 10 indicating that most of the line is emitted 
in a cocoon surrounding the pulsar, with a size of a few accretion 
radia. This cocoon is ionized and corresponds to the inner region 
of the hydrodynamical tail 

The parameters of the IGR J 17252 - 3616 are very similar 
to these of Vela X-l, excepting for the smaller wind velocity. 
We argue that the persistently large absorption column density 
is related to the hydrodynamical tail, strengthened by the low 
wind velocity. The tail is a persistent structure dissolving on a 
time scale comparable to the orbital period. 

Our interpretation can be tested using numerical hydrody- 
namical simulation and high resolutions optical/infrared spec- 
troscopy. If confirmed, it may turn out that half of the persistent 
sgHMXB have stellar wind speeds several times lower than usu- 
ally measured. 

Acknowledgements. This research has made use of NASA's Astrophysics Data 
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